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Abstract 

00 ' Selected applications of perturbative Quantumchromodynamics (QCD) to predictions of the Standard 

C^ . Model for processes at high energy colliders are reviewed with emphasis on past successes and future 

problems. This is a personal retrospective is not intended to be a comprehensive review of the field. 

I ■ 1 Introduction 

^. 



It is with great pleasure that I dedicate this contribution to Professor Piyare Lai Jain, a pioneering faculty 
member in experimental high energy and heavy ion physics at the University at Buffalo over the past fifty 
years, and a respected senior colleague of mine in the Department of Physics over the past twenty-five years. 
^^ I His mentoring, friendship and collaboration Jj have been invaluable to me. 

oo 

^ ; 2 Renormalization Prescription Dependence of QCD Predictions 



My earliest work in the field of perturbative QCD with Celmaster[2l|3] was on the renormalization prescrip- 
tion of the effective coupling as and its experimental consequences. 

Quantumchromodynamics (QCD) is an asymptotically-free quantum field theory. The scaling behavior of 
f**^ , Green's functions in the deep Euclidean region can be computed as a perturbation series in the renormalized 

coupling constant as , which depends on the characteristic energy scale involved in the process. This property 
of QCD was first used to obtain predictions for scale breaking of the naive parton model in deeply-inelastic 
lepton-hadron scattering. The domain of applicability of perturbative QCD has since been extended to 
include e+e~ annihilation to hadrons and numerous other infrared safe observables at high energy colliders[4]. 
A typical prediction of perturbative QCD is an n-th moment of a structure function in deeply-inelastic 
scattering [5] 

M„=g''"[l + fe„g2+0(/)]A„, (1) 

where g = y/^na^ is the SU(3)c coupling parameter, a„ and 6„ are calculable numbers, and the normalization 
An cannot be calculated in perturbation theory. The coefficient &„ depends on how ultraviolet divergences 
are renormalized, and so implicitly does the coupling g. The renormalization-prescription dependence of 6„ 
must be related to that of g because physical observables cannot depend on how the theory is renormalized. 
To exhibit this dependence, express g in terms of the coupling g' in a different renormalization scheme: 

g = g'[l + ag'' + 0{g'*)], (2) 
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where a is a calculable constant. To obtain the coefficient 6Jj in the second scheme, substitute the expression 
for g in Eq. Q in Eq. ((T]): 

M„ = {gT" [1 + (bn + aan)g'' + 0{g'^)]A^ = {gT" [1 + bW' + 0{g'^)]A^ . (3) 

If the perturbation series in Eqs. (jll3p were known to all orders in the couplings g, g' , then the renormalization- 
prescription dependence would have no observable consequences. In practice, the series cannot be computed 
beyond the first two or three orders of perturbation theory for most observables. The truncated series do 
depend on renormalization prescription, and differ by terms of the leading uncomputed order in the series. 
Since a^ ~ 0.1—0.3 at currently-accessible energy scales, the discrepancies are experimentally significant. 

This type of dependence on the re-definition of a small expansion parameter will arise in any truncated 
series expansion. In particular, perturbative series in the fine structure constant a of Quantumelectrody- 
namics (QED), are also subject to truncation ambiguities. The issue is less severe in QED however for two 
reasons: First, the expansion parameter a ~ 1/137.034 is much smaller than as, and so are truncation 
ambiguities. Secondly, there is a naturally preferred definition of a which is experimentally measurable in 
principle. In fact, the cross section for Compton scattering of a photon from a spin one-half charged particle 
can be shown to reduce to the classical Thompson cross section to all orders in perturbation theory in the 
linrit of zero photon energy. There is no such natural definition in QCD, cither in the low-energy limit where 
the theory becomes strongly interacting, or in the high-energy limit where the effective coupling tends to 
zero due to asymptotic freedom. 

Several approaches have been proposed to deal with the prescription dependence of predictions in per- 
turbative QCD: 

• Good universal prescriptions: The strategy is to find and use a prescription that appears to give 
a reasonably convergent perturbation series in a wide variety of processes. The archetypical example 
of this approach is the widely-used MS scheme recommended by Bardeen, Buras, Duke and MutaJ^ in 
which dimensionally-regularized quantities are renormalized by subtracting poles and certain associated 
constants in the combination 

^-^ + [log(47r) - 7e] , (4) 

where A'^ is a continuation of 4 space-time dimensions and 7e is Euler's constant. This is a modification 
of the minimal subtraction (MS) scheme of t'Hooft and Veltman[5] in which only the pole in A^ — 
4 is subtracted. Another proposal by Brodsky, Lepage and Mackenzie[7] uses light-quark vacuum 
polarization insertions to determine the scale of as in any process that does not involve gluon-gluon 
couplings in leading order. 

• Fastest apparent convergence: It is by no means guaranteed that a good universal prescription 
must exist. A different strategy suggested by Grunberg[8] tailors the prescription to each individual 
process in order to make the higher order terms in the series as small as possible. For example, in 
Eq. ([3]), one might choose 6^ = 0, so that the QCD prediction is given by the leading order term. The 
different coupling parameters, determined experimentally for example from the different moments M„, 
are related to one another by Eqs. (|2l3p . 

• The principle of minimum sensitivity: Stevenson |5] suggested choosing a renormalization pre- 
scription optimally for each process considered by requiring that the truncated perturbative expression 
be an extremum with respect to variations in the parameters of the renormaliztion group of Gell-Mann 
and Low[TD], and Stucckelberg and Petermann|llj. At next-to-leading order in QCD, predictions 
depend on one renormalization group parameter, the scale /i of the effective coupling as{fJ-). 

Each of these strategies has its advantages and shortcomings. There is no rigorous justification for any 
of them within the framework of quantum field theory. Each of them is a more or less educated guess 
concerning the size of uncomputed terms in the perturbation series. The only way to verify this guess is to 
compute these unknown terms! 



2.1 The momentum-space subtraction prescription 

Celmaster and Gonsalves[51 [3] proposed an alternative to the MS and MS prescriptions called momentum- 
space subtraction (MOM). In this approach, one considers the basic ofF-shell Green's functions that must be 
renormalized in any perturbative calculation beyond tree level. At one-loop level, these are the three-gluon, 
the quark-gluon, and (in covariant gauges) the ghost-gluon scattering amplitudes. The theory is renormalized 
by subtracting these fundamental building-block amplitudes at an energy scale that is typical of the process 
under consideration, which also serves as the scale of the effective coupling as ■ 

The motivation for the MOM prescription is the intuition that ultraviolet behavior is dominated by the 
off-shell Green's functions that need to be renormalized. If these Green's functions are primarily respon- 
sible for the divergent behavior of the unrenormalized theory, then subtracting their contributions at the 
typical energy scale of the process will minimize their effect on the perturbative correction terms. Because 
subtractions are performed on a small number of generic Green's functions, momentum-space subtraction is 
an example of a good universal prescription. 

The MOM prescription proposal suffers from some ambiguities. There are three different vertices at one- 
loop level that can be chosen to define the subtraction. Each of these off-shell vertices is gauge dependent, 
and so a particular gauge must be chosen to define the subtraction. And finally, each three-point vertex 
is has two independently- variable momenta flowing through it, and these must be fixed in relation to the 
energy scale. 

A detailed study of these vertices showed that it is possible to define a class of MOM prescriptions that 
gives reasonably-sized next-to-leading order QCD corrections in a number of processes. The corrections 
are only weakly gauge-dependent if a covariant gauge is used with the gauge parameter allowed to vary 
quite freely in the range of the usual Feynman and Landau gauges. We choose a natural and intuitively 
reasonable choice of momenta at which to subtract a three-point function, namely the symmetric Euclidean 
point Pi=P2=P3 = — M^- The off-shell Green's functions are related by Ward-Takahashi identities, which 
relate the different prescription in the MOM class to one another. For symmetric-point subtraction, the class 
of MOM prescriptions is also only weakly dependent on the chosen vertex. 

2.2 The Celmaster-Gonsalves relation for As in different schemes 

Celmaster and GonsalvesjSj showed that A's defined in different schemes can be related exactly by performing 
a one-loop calculation. 

It is well known that the renormalization group equation 

^i^ = (3{g(p)) = -(3og' - p^g" + . . . , (5) 

where fi is the QCD beta function with leading coefficient /3o = H — 2nf/3, can be integrated to express the 
scale dependence of the QCD coupling 

2, ^ ^ 1 _ /?ll0gl0g(A.VA^) 

^^^' /3olog(MVA2) Pl\og\^i^/K^) ^' 

in terms of a renormalization-group invariant QCD parameter A, which does however depend on the renor- 
malization prescription used to define g{n). This relation can be solved perturbatively for the ratio 

(7) 

where the coupling parameters are related by Eq. ^. Because A and A' are scale- invariant constants, the 
limit fi —r oo can be taken on the right hand side to give 

A^ ^ '"P [Wo\ ' ^'^ 
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where the constant a relates the two prescriptions in one-loop order. This result is exact to all orders in 
perturbation theory because of asymptotic freedom [12l [13] g{oo) = 0. 

The Celmaster-Gonsalves relation has been used to relate ALattice as measured in lattice QCD to Aj^ 
determined from high energy phenomenology [l4|, 1 1 5 j . 



3 QCD Radiative Corrections in Electron-Positron Annihilation 

On of the most powerful probes of the short-distance structure of hadrons is the virtual photon produced 
by annihilation of an energetic pair in colliding beams of electrons and positrons. This process provides a 
plethora of measurable inclusive cross sections and distributions that have been used to discover new quark 
species, and to study bound states of heavy quarks, the densities of partons in hadrons, and the fundamental 
structure of the quark-gluon interaction. The simplest observable is the total cross section for annihilation 
to hadrons, which can be expressed in the form 
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(9) 



where the sum is over quark flavors / = u,d, s,c,b,t, with Qf being the fractional quark charge, and the 
next-to-leading order coefhcient K is calculable in perturbative QCD. This simple form holds in regions far 
from flavor thresholds with the flavor sum being restricted to flavors with Ami < i?cm- The cross section is 
proportional to the number of quark colors and thus provides a direct measurement of the 3 in SU(3). The 
fractional electric charge of each new quark species is directly measured by the increase of the flavor sum as 
the threshold is crossed, and is most dramatically illustrated by the data in Fig. ([T]) on the ratio 
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Figure 1: Data on the R ratio in e+e annihilation from the PDG Review of Particle Physics [TB]. 



The leading a^/Tr correction was first computed by Schwinger|17j. The next-to-leading order coefficient 
K was computed analytically by Celmaster and Gonsalves [19l [20] . 
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where ( is Riemann's zeta function. This result agreed precisely with an earlier numerical calculation by 
Dine and SapirsteinTS] on gauge-invariant subsets of diagrams and their sum. The analytic result was also 
obtained in an independent calculation by Chetyrkin, Kataev and Tkachov[21j. 

The next-to- leading order coefficient can be expressed in other renormalizaton schemes using Eq. ^ : 

^MOM = ^ -^Cl2(f)-|n/-2/3oC(3)^ -2.193 + 0.162n/, (12) 

365 11 B 

^MS = ^ -Y^n/ + y[log(4^)-7E-4/3oC(3)]^ 7.359 -0.441n/, (13) 

where CI is Clausen's function, and the MOM scheme is defined by the quark-quark-gluon vertex in Landau 
gauge. Note that MOM and MS give a well-behaved perturbation series. The simple MS scheme is less well 
behaved. 

The theoretical prediction for R has since been extended to 0{al) by Gorishnii, Kataev and Larinp6]. 
and Surguladze and Samuel p7]. and work is in progress p5] on the 0{al) coefficient. 

3.1 The photon propagator in QCD 

The total annihilation cross section can most easily be obtained from the 2-point correlation function for the 
electromagnetic vector potential 

D^,{q) = J dxe''i-^{T[A^{x)A,{0)]) = ^ fg^, - '^\ D{-q^) + gauge terms . (14) 

The total cross section is given by the imaginary part of the photon propagator 

a(e~^e~ -^ hadrons) = —^mD(-q^) , (15) 

q 

with q"^ = -E-c.m.- 

If QCD with massless quarks were a finite theory, then the dimensionless function D{—q^) would be a pure 
number independent of q^ . Because of ultraviolet and infrared divergences, the theory must be regularized 
to compute D{—q^) as a perturbative series in the QCD couphng Us- This is most conveniently done by 
analytically continuing D{—q^) to the region of spacelike g^ < and to space-time dimension A^ = 4 — e. 
The unrenormalized function D is then real and finite, with divergences manifest as poles in e 

oo .. 

i?(-g2) = l + ^^(V)-'"rf„(.)(^)"" +Oial), (16) 

n—1 

where ao and aso are the unrenormalized (bare or zero order) QED and QCD coupling parameters respec- 
tively, and the coefficient 

, "n,n "n.n— 1 "n.l , ^^/ n /t7\ 

dn = —^ + „_i H \ hrf„,o + C(e), (17) 

is got by computing 7i-loop Feynman diagram contributions to the photon propagator. The dependence of the 
unrenormalized propagator on the sole dimensional quantity q^ in massless QCD follows from dimensional 
analysis. 

The propagator is made finite in the limit e — » by subtracting divergences at a renormalization scale 
/I. The bare QCD coupling a^o is then replaced by the effective couple as{ij). Since the electromagnetic 
current is conserved, the bare QED coupling ag can simply be replaced by its renormalized value (which 
differs from the usual fine structure constant a by terms of order a^). It is then straightforward to express 
the coefficient K in Eq. ^ in terms of the residues of the poles in Eq. p?|) . 
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Figure 2: Feynman diagrams which contribute to the photon propagator in QCD through order a^. 

3.2 The Gegenbauer expansion technique 

The residues dn^i in Eq. (|17p were computed by evaluating the Feynman diagrams in Fig. ([2]) assuming 
massless quarks and using dimensional regularization. The calculation was done by generalizing a well- 
known Chebyshev polynomial expansion technique developed by Rosner and others 22, 24 for QED to iV- 
dimensions. The appropriate generalization is provided by Gegenbauer polynomials C^ (x) with A — {N—2)/2 
and X the cosine of the polar angle in A^-dimensional spherical coordinates. 

The Gegenbauer polynomials [53] can be generated from Cq = 1 using the recursion formula 

2(n + X)xC^{x) = {n+ l)C^,+,{x) + (1 - <5„,o)(« + 2A - l)C^_i(a;) , (18) 

and obey the orthogonality relation 

' d.(i ~ x^)^c,^(x)c^(.) ^ ^».^ "j(7+ amfIaf ■ ^^^^ 

The basic idea is to expand propagator denominators of massless quarks or gluons using an TV-dimensional 
generalization of the Legendre polynomial expansion 

1 1 ^ „ /A; 



= I^E/"(t:^)^"(^i-^2)' (20) 



{ki - /C2)2 «> n=0 

where the coefficient functions /„ are related to hypergeometric functions and can be expressed as a power 
series in the dimensional regularization parameter e 
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Integrations over loop momenta can then be performed using the orthogonality relation in Eq. p9p and 
various other properties of the Gegenbauer polynomials. 



The calculation of Chetrykin, Kataev and Tkachov[21j used a similar polynomial expansion technique in 
position space [25]. 

4 The Quark Electromagnetic Form Factor in QCD 

A different method for computing the total annihilation cross section discussed in the preceding section 
would be to calculate the amplitudes for a virtual photon to decay to to 2, 3, or 4 partons, and then sum the 
individual cross sections integrated over the phase space of final state partons. Fig. ^ shows the Feynman 
diagrams which contribute to the 2-parton final states through order a^, i.e., to the quark electromagnetic 
form factor. 
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Figure 3: One-particle irreducible Feynman diagrams which contribute to the quark electromagnetic form 
factor in QCD. Solid, wavy, and dashed lines represent quark, gluon, and ghost propagators, respectively. 
Diagram IV is the only non-planar diagram in this set. 

The form factor is ultra-violet finite (when all QCD counterterms are included) because the electro- 
magnetic current is conserved, but it is infrared divergent due to soft and coUinear singularities. These 
singularities are canceled by contributions from 3- and 4-parton final states so the total cross section is 
finite. The form factor is therefore an unphysical part of a measurable physical observable. However, it is 
an object of great theoretical interest because the soft and collinear singularities are generic in that they 
determine double-logarithmic Sudakov effects in many measurable cross section — see the review article by 
Collins[29]. It was therefore of interest to study these singularities. 

A calculation [30] of all two-loop integrals required for an analytical calculation of the form factor was 
completed in 1983, including all terms of ©(e"") with < ti < 4. The most challenging sets of integrals to 
compute were those which contribute to the box-plus-triangle diagram III and the non-planar diagram IV 
in Fig. ([3]). The calculations were done by converting the two-loop momentum integrals to 5-dimensional 



integrals over Feynman parameters. The Feynman parameter integration were performed analytically in suc- 
cession using various techniques including infinite-series expansion of binomial expressions where necessary. 
Care had to be taken to ensure that the binomial expansions converged after each successive integration: 
where necessary, the integration interval was subdivided and appropriately converging expansions used in 
each region. 

The results for the two-loop scalar integrals were verified by van Neerven[Jl] and by Kramer and 
Lampe[32]. Gehrmann, Huber and Maitre[33] have recently obtained elegant exact expressions for the 
integrals in diagram IV in N dimensions. 

These two- loop results have been used in several 0{al) QCD calculations including the Drell-Yan total 
cross section by Hamberg, van Neerven and Matsuura[31], the total cross section for Higgs production by 
Harlander and Kilgore[3S] and Anastasiou and Melnikov[36j. the Drell-Yan rapidity distribution by Anasta- 
siou, Dixon, Melnikov and Petriello^S^, and the Higgs total cross section and Drell-Yan and Higgs rapidity 
distributions by Ravindran, Smith and van Neerven [38]. 

In a closely related calculation f395, a tensor decomposition scheme that is exact in N dimensions was 
used to compute the one-loop QCD corrections to the e'^e" total annihilation cross section. Unlike the 
process 7* —^ qq in which there is a single scalar form factor, the process 7* —^ qqg involves 6 independent 
invariant tensors and corresponding form factors in A'^ = 4 dimensions. When the amplitude is continued to A^ 
dimensions, there is a set of 17 linearly-independent tensors and corresponding form factors into which it can 
conveniently be decomposed. The decomposition makes it possible to compute the dimensionally-regularized 
amplitude very efficiently in ©(a^), and the technique readily extends to two or more loops. 

Over the past few years, a number of results on form factors have appeared in the literature, including 
analytic results for arbitrary masses and momenta in QEDjJDlST] and in QCD [151133], and impressive new 
results on massless quark and gluon form factors at three loops '44]. 

5 Chromo-Electroweak Interference and Parity Violation 

In the Standard Model the process qq -^ qqV (where V = W^, Z^ or 7) can occur via gluon exchange and 
also via W or Z^ exchange. The corresponding chromodynamic and electroweak amplitudes can interfere 
with one another [45] as shown in Fig. (|4]). These interference cross sections are largest when the exchanged 
W^ or Z^ is on-shell when they are also odd under parity. Interference cross sections computed using 
helicity- amplitude techniques were analyzed [3^ for all interesting subprocesses as well as for the processes 
qq — > qqll in which the lepton pair II comes from the decay of V on-shell. Parity-violating asymmetries were 
defined and presented at the parton level and for the hadronic processes pp or pp ^ V + 2 jets or Z^ -I- 2 
jets, see Fig. ([S]). These asymmetries are independent of the polarizations of all particles involved, and do 
not require that the flavors of the jet partons be measured. They are generally of order 0.01 pb at energies 
Vs > 1 TeV. 
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Figure 4: A solid line represents a quark or antiquark with color c and flavor /. The wavy line represents a 
color-neutral electroweak boson and the curly line a flavor-neutral gluon. The two amplitudes can interfere 
because their initial states (on the left) can have identical quantum numbers as can their final states (on the 
right). 

A simple way of defining experimental observables that are sensitive to this parity-odd character of the 
interference terms is as follows: Imagine that the incident beams lie in the plane of a mirror, as in Fig. ^. 
If the incident beams are not polarized, the initial state is invariant under reflection in this mirror. We will 




Figure 5: Hadrons hi and /12 collide to produce an on-shell electroweak boson V and two hadron jets ji 
and j2 at large transverse momentum. 

also assume that the spins of the final state particles are not detected, i.e., particles are identified by their 
momenta and internal quantum numbers only. A parity-odd contribution to the cross section can make the 
probabilities for observing a particular event (i.e., a particular configuration of final state particles) and its 
mirror image different from one another. Since the events are continuously distributed in phase space, the 
likelihood of finding an event and its geometrical mirror image in any finite sample of events is vanishingly 
small. To decide experimentally whether or not there is an asymmetry with respect to mirror reflection in 
the event sample, one must count the number of events that fall in some region of phase space (which we 
will call a "bin" ) and the number of events that fall in the mirror image of this region (which we will the 
"image bin" ) and then decide whether or not there is a statistically significant difference between these two 
numbers. This difference can be compared with a "parity-violating asymmetry" which we define as follows: 



aP^(bin) ^ da 

bin 



dcT 



(22) 



image bin 



The integral over the bin includes an implicit sum over all final state quantum numbers that are not ob- 
served, i.e., color quantum numbers and any spin or flavor quantum numbers that cannot be experimentally 
measured. 

A set of bins will have to be judiciously chosen to maximize the observed effects of parity violation, i.e., 
to yield the largest cumulative asymmetry. Thus we also define a cumulative parity-violating asymmetry as 
follows: 

^Pv ^ ^ I „pv(bin) I (23) 

bins 

where the sum is taken over bins that do not overlap with one another. A theoretical upper bound on Ap^ 
is given by 
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These asymmetries might be observable in pp and pp collisions above the threshold for production of pairs 
of electroweak bosons, i.e., at center of mass energies in the TeV range and these asymmetries are comparable 
in magnitude to the pair-production cross sections. Fig. ([7]) compares the parity-violating signal with the 
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Figure 6: (a) Momenta in the x—z plane, pi is the momentum of the quark (parton-level asymmetry), 
the proton (pp colhsions), or one of the protons chosen arbitrarily (pp collisions), ps is the momentum 
oiV = W^,Z^ (3-particle final state) or of the lepton (4-particle final state). A right-handed coordinate 
system is defined such that pf > and p| — 0. (b) Momenta in the x—y plane, p^ and p4 are the momenta 
of the two jet partons. Parity is violated if the event shown on the left and its mirror image shown on the 
right occur with different probabilities. 

QCD background [5D] in W^ + 2 jet production. Unlike the W^-pair production cross section which is parity 
conserving, the interference contribution is parity-violating, and this might make it easier to observe above a 
rather formidable QCD background of electroweak boson -|- 2 jet events. Actual observation of these effects 
will require somewhat higher integrated luminosities than are currently available for example at the Tevatron 
at Fermilab. It remains to be seen whether these subtle Standard Model predictions will be observable at a 
high-energy multi-TeV hadron collider such as the LHC. Analogous parity-violating signatures may also be 
observable in electron-positron annihilations [IJ at the ILC. 



6 Electroweak Boson Production at Large Transverse Momentum 

The production of electroweak bosons at large transverse momentum is one of the most important processes 
at current and future hadron colliders. The QCD-improved parton model predicts that if the intrinsic energy 
scale involved is sufficiently large, the inclusive cross section for the process 



hi{Pi) + h2iP2)-^V{Q)+X 



(26) 



where hi, i = 1,2 are unpolarized hadrons with momenta Pi, can be reliably computed using the following 
approximate factorized form: 
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Here Eq = Q°, a and h stand for quarks, antiquarks or gluons, f'^{x,M'^) is the probability density for 
finding parton a with momentum fraction x in hadron h if it is probed at scale M"^, and a°'''{pi,p2, M'^) is 
the perturbative cross section for the process 



a{pi) + b{p2)^V{Q)+X , 



(28) 
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Figure 7: Comparison of signal to background in W~ + 2 jet production, (a) Binned cross sections with 
200 GeV < ^/Sl2 < 600 GeV. The parity-violating asymmetries have been multiplied by the factors (100 
and 1000) indicated to show them on the same scale as the QCD background, (b) Variation of the signal to 
background ratio in pp and pp collisions as a function of the colliding beam energy. 

from which coUinear singularities arising from radiation off massless partons have been factorized out at 
scale M^ and impHcitly included in the scale-dependent parton densities /^ (a;,M^). 

A complete analytical calculation of the next-to-leading order QCD radiative corrections to the inclusive 
cross sections parton + parton —> V + X, where V is an on-shell W or Z^ with transverse momentum Qt 
of order Mw, or a massive virtual photon with Qt of order of its invariant mass, was done by Gonsalves, 
Pawlowski and Wai[49 and by Arnold and Reno 50 . This work completed an earlier calculation by Ellis, 
Martinelli and PetronzioL48j of the non-singlet (NS) contributions to lepton-pair production at large Qt- 
Numerical predictions for W, Z and 7* production at collider energies were computed, and the dependence 
of the radiative corrections on the choice of renormalization and factorization scales was studied. Results [HT| 
on varying the renormalization and factorization scales independently were used to compare the FAC[5] and 
PMS|n] prescriptions. These results showed that the QCD-improved parton model can be used to make firm 
and reliable predictions for electroweak boson production at large Qt- 

In Fig. ([5]), the theoretical predictions are compared with measurements [5^ [55] of the W transverse 
momentum distributions by the UAl and UA2 collaborations at the CERN SppS collider. It is clear that 
data are consistent with the QCD predictions, but that they are not accurate enough to discriminate between 
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Figure 8: W production in pp collisions at 630 GeV. The solid, dashed, and dotted lines represent the Full 
(MRSB), Full (MRSE), and NS (MRSB) predictions. The corresponding total cross sections a are 5.8, 5.1, 
and 3.5 nb, respectively. 

the MRSB and MRSE distributions [54] or between the singlet and non-singlet theoretical predictions. There 
is a tantalizing, but statistically insignificant, hint in the data that the experimental cross section might be 
larger than the theoretical prediction at large values of Qt- 

Fig. ([9|) shows a comparison of theoretical predictions with experimental measurements by the DO Col- 
laboration. The theoretical results in this figure have been improved at small Qt by resumming Sudakov 
logarithms of Q^/M^^. The data are completely consistent with the predictions and there is no evidence in 
the tail of the distribution for new physics beyond the Standard Model at large Qt- 

6.1 Large logarithmic soft-gluon corrections 

Predictions for the W transverse momentum distribution at next-to-leading order (NLO) in a^ and including 
the next-to-next-to-leading order (NNLO) soft-gluon corrections at LHC energies were presented recently by 
Gonsalves, Kidonakis and Sabio-Vera[56j. 

Let S2 be the invariant mass of the final state partons recoiling against the W boson. In general, the 
partonic cross section a includes distributions with respect to S2 at n-th order in the QCD coupling ag of 
the type 



ln"(52/Q|) 



S2 



defined by their integral with any smooth function / by 



/" 

Jo 



dS2/(s2) 



In™(s2/Q|^) 



S2 



m < 2n — 1 



as2'^L:i^^l^[fis2)-fm 



(29) 



m 



1 



In 



S2 
rn + l I ^ 2 max 
■^T 



/(o). 



(30) 



These "plus" distributions are the remnants of cancellations between real and virtual contributions to the 
cross section. Let us make use of the terminology that at n-th order in a^ the leading logarithms (LL) 
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Figure 9: W production in pp collisions at 1.8 TcV measured by the DO CoUaborationf 



are those with m = 2n — 1 in Eq. ((29| . the next-to-leading logarithms (NLL) are those with m = 2n — 2, 
the next-to-next-to-leading logarithms (NNLL) are those with m — 2n — 3, and the next-to-next-to-next- 
to-leading logarithms (NNNLL) are those with m = 2n - 4. The symbol "NNLO-NNNLL" means that 
soft-gluon contributions through NNNLL to the NNLO corrections have been included and added to the 
complete NLO result. 

Fig. pni shows the differential cross section da/dQ^ at high Qt with ^/S = 14 TcV for the two scale 
values, QtI'2 and 2Qt^ often used to display the uncertainty due to scale variation. Note that while the 
variation of the Born cross section is significant and the variation at NLO is similar to LO, at NNLO-NNNLL 
it is very small. In fact the two NNLO-NNNLL curves lie very close to or on top of each other. With this 
calculation, reliable predictions of perturbative QCD, including soft-gluon effects, for the W transverse 
momentum distribution are available and await comparison with experimental measurements at the LHC 
in the near future. Any deviation from these predictions would be an unambiguous signal of new physics 
beyond the Standard Model. 

It is obviously worthwhile to compute the complete NNLO corrections to the transverse momentum 
distribution. This would further reduce the scale dependence of the predictions, and it would also provide 
information on soft and logarithms at 0{a\). Fixed order QCD calculations for analogous processes in e+e" 
annihilation have been completed over the past few years, see e.g. Gehrmann-De Ridder et al.j57j, and 
references therein. An ©(a^) calculation at finite Qt would require the two-loop corrections to subprocesses 
such as qq — > W^5, one-loop corrections to subprocesses such as qq — > Wgg, and tree-level amplitudes such 
as qq -^ Wggg. 
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Figure 10: The differential cross section, da/dQ-^, for W production in pp collisions at the LHC with ^/S = 14 
TeV and fi = ^ip = f^R = Qt/2 or 2Qt. Shown are the LO, NLO, and NNLO-NNNLL results using the 
MRST 2002 LO, NLO and NNLO parton distributions 64J . The upper hues are with fi = Qt/^, the lower 
lines with fi = 2Qt- 

7 Heavy Quark Triangle Diagram Contributions at High Energies 

Because the masses of the W and Z bosons are large compared with the masses of the quark and gluon 
constituents in the colliding hadrons, it is generally a good approximation to assume that the quarks are 
massless in computing the production cross sections for these electroweak bosons at high energies. An 
exception, of course is the top quark, whose mass mt = 174 GeV, is approximately twice as large as that of 
the W or Z. 

The top quark can contribute to production amplitudes in three ways: as an initial state partonic 
constituent of one of the colliding hadrons, as an on-shell particle in the final state, and as a virtual particle. 
The density of top quarks in the colliding hadrons can safely be neglected. The production of real top quarks 
can safely be treated in a threshold approximation with the mass taken to be infinite below threshold and 
zero above. Top quarks in loops pose a delicate problem because the Standard Model is renormalizable only 
if the left-handed quarks occur in 811(2)^ doublets. This requirement is imposed by the well-known triangle 
anomaly (which was discovered in QED by Adler[55j and Bell and Jackiw[55]) in the effective coupling of the 
axial- vector current to two gluons: if the top quark is omitted from the loop the amplitude is divergent in a 
non-renormalizable way; if the top quark is included and assumed massless the top and bottom contributions 
cancel; and if the top quark is included with a finite mass the contribution grows as the logarithm of the 
ratio of top to bottom masses. 

The contributions in second order QCD of the diagrams in the left panel of Fig. (fTTj) involving a 
heavy-quark triangle with vector, vector and axial-vector vertices were studied by Gonsalves, Hung and 
Pawlowski[62]. The contributions from the process qq — *■ Zgg ((a) in Fig. pT|) ) had been analyzed earlier by 
Dicus and Willenbrock [61] . who concluded that the contributions were numerically small at Tevatron ener- 
gies. The right panel of Fig. (|lip shows the numerical magnitudes of the various subprocess contributions 
at Tevatron energies. 

Rijken and van Neerven[53] subsequently computed all heavy quark contributions to the W and Z 
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Figure 11: Left panel: Diagrams which contribute to the total cross section for producing a Z boson in 
hadron collisions. Q and q represent heavy and massless quarks, respectively, and g a gluon. Propagators 
cut by a dotted line are taken to be on shell in the initial or final state. Right panel: Contributions from 
heavy-quark triangle diagrams to the total cross section for Z production in pp collisions at 2 TeV. m is the 
mass of the heavy quark, which is assumed to have weak isospin +i. A quark with weak isospin — i would 
contribute with opposite sign. The total cross section at 2 TeV is approximately 5 nb. 



total cross sections, neglecting only diagrams with heavy quarks in the initial state. They found that the 
contributions from the top quark were negligible at Tevatron energies, but would contribute significantly at 
LHC energies with magnitude comparable to the 0[a\^ QCD corrections. 

It is important to evaluate the top-quark-mass contributions to other measurable cross sections at LHC 
energies, in particular the W and Z transverse momentum distributions. Fig. (fT^ shows the Z transverse- 
momentum distribution at Tevatron and LHC energies including quark-triangle diagrams shown in Fig. (jlip . 
Quark masses are taken into account using a simple threshold prescription ^49j by inserting a step function 
^(Qt ~ 4mp into sums over quark flavors / in the subprocess cross sections. There is a significant change 
in the shape of the distribution at LHC energies due specifically to the virtual triangle- anomaly diagrams. 

Because the triangle-anomaly contributions are comparable in magnitude to the 0(c?^^ QCD radiative 
corrections, it is essential to include the top and bottom quark masses in the virtual triangle diagrams. A 
more precise calculation of the triangle-anomaly contributions to the Z transverse momentum distribution 
at LHC energies will be published elsewhere. 
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Figure 12: Z production at the LHC and the Tevatron. The data points represent NLO predictions 
computed using MRST 2002 NLO parton distributions 64 . Note the change in shape of the distribution at 
LHC energies (upper curves) compared with Tevatron energies (lower curves) when the anomaly diagrams 
are included. The anomalous contributions are negligibly small at the Tevatron. 

References 

[1] p. L. Jain, G. Singh and R. Gonsalves, Mod. Phys. Lett. A 12 (1997) 1101. 

[2] W. Celmaster and R. J. Gonsalves, Phys. Rev. Lett. 42 (1979) 1435. 

[3] W. Celmaster and R. J. Gonsalves, Phys. Rev. D 20 (1979) 1420. 

[4] R. K. Elhs, W. J. Stirling and B. R. Webber, QCD and Collider Physics (Cambridge University Press, 
Cambridge, 1996). 

[5] W. A. Bardeen, A. J. Buras, D. W. Duke and T. Muta, Phys. Rev. D 18 (1978) 3998. 

[6] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 44, 189 (1972). 

[7] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28 (1983) 228. 

[8] G. Grunberg, Phys. Lett. B 95 (1980) 70; Erratum ibtd. B 110 (1982) 501. 

[9] P. M. Stevenson, Phys. Rev. D 23 (1981) 2916. 

[10] M. Gell-Mann and F. E. Low, Phys. Rev. 95 (1954) 1300. 

[11] E. C. G. Stueckelberg and A. Petermann, Helv. Phys. Acta 26 (1953) 499. 

[12] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30 (1973) 1343. 

[13] H. D. Politzer, Phys. Rev. Lett. 30 (1973) 1346. 

[14] A. Hasenfratz and P. Hasenfratz, Phys. Lett. B 93 (1980) 165. 

[15] R. F. Dashen and D. J. Gross, Phys. Rev. D 23 (1981) 2340. 



16 



[16] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[17] J. S. Schwinger, Phys. Rev. 74 (1948) 1439. 

[18] M. Dine and J. R. Sapirstein, Phys. Rev. Lett. 43 (1979) 668. 

[19] W. Celmaster and R. J. Gonsalves, Phys. Rev. Lett. 44 (1980) 560. 

[20] W. Celmaster and R. J. Gonsalves, Phys. Rev. D 21 (1980) 3112. 

[21] K. G. Chetyrkin, A. L. Kataev and F. V. Tkachov, Phys. Lett. B 85 (1979) 277. 

[22] E. De Rafael and J. L. Rosner, Annals Phys. 82 (1974) 369. 

[23] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover, New York, 1972) 

[24] M. J. Levine and R. Roskies, Phys. Rev. D 9 (1974) 421. 

[25] A. A. Vladimirov, Theor. Math. Phys. 43 (1980) 417 [Teor. Mat. Fiz. 43 (1980) 210]. 

[26] S. G. Gorishnii, A. L. Kataev and S. A. Larin, Phys. Lett. B 212 (1988) 238. 

[27] L. R. Surguladze and M. A. Samuel, Phys. Rev. Lett. 66 (1991) 560; Erratum ihid. 66 (1991) 2416. 

[28] P. A. Baikov, K. G Chetyrkin and J. H. Kuhn, Nucl. Phys. Proc. Suppl. 157 (2006) 27 
|arXiv:hep-ph/060212^ . 



[29] J. C. Collins, Adv. Ser. Direct. High Energy Phys. 5 (1989) 573 |arXiv:hep-ph/0312336| . 

[30] R. J. Gonsalves, Phys. Rev. D 28 (1983) 1542. 

[31] W. L. van Neerven, Nucl. Phys. B 268 (1986) 453. 

[32] G. Kramer and B. Lampe, J. Math. Phys. 28 (1987) 945. 

[33] T. Gehrmann, T. Huber and D. Maitre, Phys. Lett. B 622 (2005) 295 | arXiv:hep-ph/050706 11 . 

[34] R. Hamberg, W. L. van Neerven and T. Matsuura, Nucl. Phys. B 359 (1991) 343 Erratum, ihid. B 644 
(2002) 403. 

[35] R. V. Harlander and W. B. Kilgore, Phys. Rev. Lett. 88 (2002) 201801 |arXiv:hep-ph/0201206| . 

[36] C. Anastasiou and K. Melnikov, Nucl. Phys. B 646, 220 (2002) [arXiv:hep-ph/0207004| . 

[37] C. Anastasiou, L. J. Dixon, K. Melnikov and F. Petriello, Phys. Rev. D 69 (2004) 094008 
|arXiv:hep-ph/0312266 . 

[38] V. Ravindran, J. Smith and W. L. van Neerven, Nucl. P hys. B 665, 325 (2003) |arXiv:hep-ph/0302135| ; 



Nucl. Phys. B 767, 100 (2007) arXiv:hep-ph/0608308 



[39] R. J. Gonsalves, Phys. Rev. D 34 (1986) 1316. 

[40] P. Mastroha and E. Remiddi, Nucl. Phys. B 664 (2003) 341 |arXiv:hep-ph/0302162| . 



[41] R. Bonciani, P. Mastroha and E. Remiddi, Nucl. Phys. B 676 (2004) 399 [arXiv:hep-ph70307295 . 

[42] W. Bernreuther, R. Bonciani, T. Gehrmann, R. Heinesch, T. Leineweber, P. Mastrolia and 
E . Remiddi, Phys. Rev . Lett. 95 (2005) 261802 [arXiv:hep-ph70509341]; Phys. Re v. D 72 (2005) 096002 
|arXiv:hep-ph/0508254; Nucl. Phys. B 712 (2005) 229 f arXiv:hep-ph/0412259 ; Nucl. Phys. B 706 
(2005) 245 [arXiv:hep-ph/0406046] . 

17 



[43] W. Bernreuther, R. Bonciani, T. Gehrmann, R. Hcincsch, T. Leineweber and E. Remiddi, Nucl. Phys. 
B 723 (2005) 91 farXiv:hep-ph70504190 . 



[44] S. Moch, J. A. M. Vermaseren and A. Vogt, Phys. Lett. B 625 (2005) 245 |arXiv:hep-ph/0508055| ; 
JHEP 0508 (2005) 049 [arXiv:hep-ph/0507039] . 

[45] R. J. Gonsalves, Phys. Rev. Lett. 56 (1986) 1647. 

[46] R. J. Gonsalves and C. F. Wai, Phys. Rev. D 49 (1994) 190; Erratum, ibid. D 51 (1995) 1428. 
I arXiv:hep-ph/961027 8 . 

[47] S. Raina, "Interference phenomena in electron positron annihilation," Ph.D. Thesis, University at Buf- 
falo, 2002, UML30-52534. 

[48] R. K. Ellis, G. MartincUi and R. Pctronzio, Nucl. Phys. B 211 (1983) 106. 

[49] R. J. Gonsalves, J. Pawlowski and C. F. Wai, Phys. Rev. D 40 (1989) 2245. 

[50] P. B. Arnold and M. H. Reno, Nucl. Phys. B 319 (1989) 37; Erratum, tbid. B 330 (1990) 284. 

[51] R. J. Gonsalves, J. Pawlowski and C. F. Wai, Phys. Lett. B 252 (1990) 663. 

[52] C. Albajar et al. [UAl Collaboration], Phys. Lett. B 193 (1987) 389. 

[53] R. Ansari et al. [UA2 Collaboration], Phys. Lett. B 215 (1988) 175. 

[54] A. D. Martin, R. G. Roberts and W. J. Stirling, Mod. Phys. Lett. A 4 (1989) 1135. 



[55] B. Abbott et al. [DO Collaboration], Phys. Lett. B 513 (2001) 292 |arXiv:hep-ex/0010026] . 

[56] R. J. G onsalves, N. Kidonakis and A. S. Vera, Phys. Rev. Lett. 95 (2005) 222001 
|arXiv:hep-ph/0507317] . 

[57] A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover and G. Heinrich, Jet rates in electron-positron 
annihilation at 0{a^g) in QCD, arXiv:0802.0813. [hep-ph]. 

[58] S. L. Adler, Phys. Rev. 177 (1969) 2426. 

[59] J. S. BeU and R. Jackiw, Nuovo Cim. A 60 (1969) 47. 

[60] R. K. Ellis and R. J. Gonsalves, in Proceedings of the Oregon Workshop on Super High Energy Physics, 
ed. D.E. Soper (World Scientific, Singapore, 1986), pp. 287-293. 

[61] D. A. Dicus and S. S. D. Willenbrock, Phys. Rev. D 34 (1986) 148. 

[62] R. J. Gonsalves, C. M. Hung and J. Pawlowski, Phys. Rev. D 46 (1992) 4930. 

[63] P. J. Rijken and W. L. van Neerven, Phys. Rev. D 52 (1995) 149 arXiv:hep-ph/9501373| . 

[64] A. D. Martin, R. G. Roberts, W. J. Stirhng and R. S. Thorne, Eur. Phys. J. C 28 (2003) 455 
[arXiv:hep-ph/0211080, . 



18 



